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Abstract 

It is shown that in 2 + 1 dimensions, a constant magnetic field is a strong catalyst of 
dynamical flavor symmetry breaking, leading to generating a fermion dynamical mass even 
at the weakest attractive interaction between fermions. The essence of this effect is that 
in a magnetic field, in 2 + 1 dimensions, the dynamics of fermion pairing is essentially one- 
dimensional. The effect is illustrated in the Nambu-Jona-Lasinio model in a magnetic field. 
The low-energy effective action in this model is derived and the thermodynamic properties of 
the model are considered. The relevance of this effect for planar condensed matter systems 
and for 3 + 1 dimensional theories at high temperature is pointed out. 



1 Introduction 

Recently there has been considerable interest in relativistic field models in 2 + 1 dimensions. Besides 
that the sophisticated dynamics of these models is interesting in itself, the models also serve as 
effective theories for the description of long wavelength excitations in planar condensed matter 
systems Jl|, |2|. Also, their dynamics imitates the dynamics of 3 + 1 dimensional theories at high 
temperature. 

In this paper, we will show that a constant magnetic field acts as a strong catalyst of dynamical 
flavor symmetry breaking (generating fermion masses) in 2 + 1 dimensions. We will in particular 
show that there is a striking similarity between the role of the magnetic field in (2 + l)-dimensional 
models and the role of the Fermi surface in the Bardeen-Cooper-Schrieffer (BCS) theory of su- 
perconductivity ||: both of them enhance the interactions of fermions in the infrared region (at 
small momenta) thus leading to generating a fermion dynamical mass (energy gap in the fermion 
spectrum) even at the weakest attractive interaction between fermions. 

We note that necessity of a supercritical dynamics (with an effective coupling constant g being 
larger than a critical value g c > 0) for generating fermion dynamical masses is a common feature of 
the dynamics in 3 + 1 and 2 + 1 dimensions ||. As will be shown in this paper, in 2 + 1 dimensions, a 
magnetic field reduces the value of the critical coupling to zero. We note that the fact that a constant 
magnetic field enhances fermion dynamical masses in the Nambu-Jona-Lasinio (NJL) model H has 
been already pointed out in Ref. ||. However, what we will show is not just that the magnetic 
field enhances the dynamical mass created by the strong (supercritical) NJL interaction but that 
in 2 + 1 dimensions, it catalyses generating the mass even at the weakest attractive interaction. 
The essence of this effect is that in a magnetic field, in 2 + 1 dimensions, the dynamics of fermion 
pairing (relating essentially to fermions at the lowest Landau level) is one- dimensional (see Sec. 2). 

We stress that this effect is universal, i.e., model independent, in 2 + 1 dimensions. This 
point may be important in connection with consideration of this effect in such condensed matter 
phenomena as the quantum Hall effect |TJ and high temperature superconductivity ||. Another, 
potentially interesting, application for this effect may be in 3 + 1 dimensional theories at high 
temperature (quark-gluon plasma in a magnetic field, for example): indeed, at high temperature, 
their dynamics effectively reduces to the dynamics of 2 + 1 dimensional theories. 



As a soluble example, we shall consider the NJL model in a magnetic field, in the leading order 
in 1/N expansion. We shall derive the low-energy effective action in the model and also study its 
thermodynamic properties. 

The paper is organized as follows. In Sees. 2 and 3 we consider the problem of a relativistic 
fermion in a magnetic field in 2 + 1 dimensions. We show that the roots of the fact that a magnetic 
field is a strong catalyst of dynamical flavor symmetry breaking in 2 + 1 dimensions are actually in 
this problem. In Sees. 4-7 we study the NJL model in a magnetic field in 2 + 1 dimensions. We derive 
the low energy effective action and determine the spectrum of long wavelength collective excitations 
in this model. In Sec. 8 we study the thermodynamic properties of the NJL model in a magnetic 
field. We show that there is a symmetry restoring phase transition at high temperature. In Sec. 9 
we summarize the main results of the paper. In Appendices A and B, some useful formulas and 
relations are derived. In Appendix C the reliability of the 1/N expansion in this model is discussed. 

2 Dynamical Flavor Symmetry Breaking in the Problem 
of Fermions in a Constant Magnetic Field 

In this section we will discuss the problem of relativistic fermions in a magnetic field in 2 + 1 
dimensions. We will show that the roots of the fact that a magnetic field is a strong catalyst of 
flavor symmetry breaking are actually in this dynamics, which plays here the role similar to that of 
the dynamics of the ideal Bose gas for an almost ideal Bose gas in the theory of superfluidity [[| . 

The Lagrangian density in the problem of a relativistic fermion in a constant magnetic field B 
takes the following form in 2 + 1 dimensions: 

' r *,(zfD^-m)*l,// = 0,l,2 



where the covariant derivative is 

D li = d ll -ieAf, Af = -Bx 2 5^. (2) 

In 2 + 1 dimensions, there are two inequivalent representations of the Dirac algebra: 

7° = cr 3 ,7 1 =ioi,7 2 =ia 2 (3) 



and 



7° 



-^3,7 X 



-id, 7 2 = -ia 2 



(4) 



where <Tj are Pauli matrices. 

Let us begin by considering the representation (|3]). The energy spectrum in the problem (|1|) 
depends on the sign of eB; let us first assume that eB > 0. Then, the energy spectrum takes the 
form (to be concrete, we assume that m > 0) [[7): 



E = uj = m, 

E n = ±u n = ±\/m 2 + 2\eB\n, n = 1, 2, 



(the Landau levels). 

The general solution is 
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Here tu n (0 = (7r 1 / 2 2 n n!)~ 1 / 2 e~^ ^ 2 H n (^), H n (£) are Hermite polynomials,/ = |e_B|~ 1//2 is the mag- 
netic length, k = 2iip/Li (p = 0, ±1, ±2, . . .), L\ is the size in the Xi-direction, £ = X2/I + fcZ. As 
Li — > 00, the density of the states at each level n is \eB\/2ir 0. 

Thus the lowest Landau level with n = is special: while at n > 1, there are solutions corre- 
sponding to both fermion (£?„ = cj n ) and antifermion (E7 n = — u n ) states, the solution with n = 
describes only fermion states. 



As eB — > —eB < 0, the solution becomes 

V( X ) = J2 Onp<-p(x) + J2 K P <p(x) (8) 

n,p np 

where the charge conjugate spinors v c and u c are v c = 72V T , u c = 72« T . Therefore at eB < 0, the 
lowest Landau level with n = describes antifermion states. 

If we used the representation (P for Dirac's matrices, the general solution would be given by 
Eq. (|5p with u np (x), v np (x) being substituted by (— l) n v n _ p (— x), (— l) n M„_ p (— x) (the factor (— l) n 
is introduced here for convenience): 



*(*) =Y,Cnp(-l) n V n -p(-x)+J2dip(-l) n U n p(-x). (9) 

n,p n,p 

We note that the mass term in the Lagrangian density (jl|) violates parity defined by 

P : #(z°, x 1 , x 2 ) -► ai*(x°, -x\ x 2 ). (10) 

However if one uses the four-component fermions || , connected with a four- dimensional (reducible) 
representation of Dirac's matrices 

V=(* ° l-^h ° )„•=(* ° i, (n; 

V -a 3 ) \ -ztn / \ -zcr 2 

the mass term in the Lagrangian density 

1 



c 2 



*, (i7 M -D M - m) * (12) 



preserves parity defined now as 



P : #(x°, x\ x 2 ) -> -7V^(^ , -x 1 , x 2 ), (13) 

z 



where the Dirac matrix 7 3 is 

o ^ 



7 3 = z 



(14) 



/ 



/ 

The important point is that the Lagrangian density ([TJ) with m = is invariant under the U(2) 
(flavor) transformations with the generators 

T = I, T x = 75, T 2 = i 7 3 , T 3 = 7 3 7 5 , (15) 

i 



where 
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i7°7 1 7 2 7 3 




(16) 



The mass term breaks this symmetry down to the U(l) x U(l) with the generators T and T 3 . 

We note that the four-component fermions appear in low-energy effective actions describing 
planar condensed matter systems with two sublattices 0. Actually, usually they appear in the 
actions without the mass term, and the important problem is to establish a criterion of dynamical 
flavor symmetry breaking which may occur as a result of interaction between fermions [4,9-16]. As 
was already indicated in Section 1, dynamical flavor symmetry breaking in 2 + 1 dimensions usually 
takes place only at a rather strong effective coupling between fermions. 

Let us now show that at m = and B / 0, the dynamical breakdown of the U(2) flavor 
symmetry takes place already in the theory (p~2"D, even without any additional interaction between 
fermions. In order to prove this, we will show that in the limit m — » 0, the flavor condensate 
(0|**|0) is nonzero: (0|M|0) = -\eB\/2it. 

The condensate (0|\l/\I/|0) is expressed through the fermion propagator S(x,y) = 
= (0\T*(x)*{y)\Q) : 

(0|M|0) = - lim trS{x, y). (17) 

x >y 

The propagator S is calculated (following the Schwinger (proper time) approach |17|) in Appendix A. 
It is 



S{x,y) 
S(x) 



x 



exp^e / Af l dz A )S{x - y) 
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(19) 



where C v = g^ + ((F art ) 2 )' fl/ (1 - eBscot(eBs))/B 2 , F*f = d^A^ - d v A ext with Aj? given in 
Eq. (0). The integral in Eq. ( |IBD is calculated along the straight line. 
The Fourier transform S(k) = J d 3 xe lkx S(x) is: 



S(k) 



dsexp 



, , 2 tan(eBs) 



—ism + isk n — isk z 



"-- eBs 
(k + m+ (&V - &V) tan(eSs)) (l + 7V tan(eBs) 
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(20) 



Transferring this expression into Euclidean space (k° —* ik^, s — > —is), we find: 



^fi(fc) — —i I dsexp 
Jo 



2 ,2 ,2 tanll ( e5s ) 



-s m^ + A;, + fe' 



e^s 



(21) 



( 73 = -17°, 7! = 7 1 ,72 



-fe|*TAi + m + -(hli - hl2) tanh(ei?s) j f 1 + -7172 ta,nh(eBs] 



7 2 are antihermitian matrices). 



From Eqs. (]17|), (0) and (|2lD we find the following expression for the condensate: 
(0|M|0) = 



: tr I d 3 kS E (k) 
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(22) 



'AV 2tt 

where A is an ultraviolet cutoff. 

Thus in a constant magnetic field, spontaneous breakdown of the flavor U{2) symmetry takes 
place even though fermions do not acquire mass (m = 0)! Note that in 3 + 1 dimensions, the 
result would be (Ol^^lO) ~ mlnm — > as m — > 0. Therefore this is a specific 2 + 1 dimensional 
phenomenon. 

What is the physical basis of this phenomenon? In order to answer this question, we note that 
a singular, 1/m, behavior of the integral in Eq. Q2"2"D is formed at large, s — > 00, distances (s is the 
proper time coordinate). Actually one can see from Eq. fl22|) that the magnetic field effectively 
removes the two space dimensions in the infrared region thus reducing the dynamics to a one- 
dimensional dynamics which has much more severe infrared singularities. From this viewpoint, the 
action of the magnetic field in the present problem is similar to that of the Fermi surface in the 
BCS theory @. 

This point is intimately connected with the form of the energy spectrum of fermions in a constant 
magnetic field. Eqs. (f|), ( [TTj) and (^) imply that for the four- component fermions, the energy 
spectrum is: 



En 



±c<Jn = ±m, 



E,, 



±u! n = ±\/m 2 + 2\eB\n, n > 1. 



(23) 



The density of the states with the energy ±ujq = ±m is \eB\/2ir, and it is \eB\j-n at n > 1. As 
m — > 0, the energy _E goes to zero and therefore there is the infinite vacuum degeneracy in this case. 
The value of the condensate ( ^2|) is equal to the density of the states at the lowest Landau level. 
This implies that spontaneous flavor symmetry breaking is intimately connected with the dynamics 
of fermions at this level. In particular, since this dynamics, described by one continuous variable 
&3 = — ik°, is one-dimensional, we get a simple explanation of the one-dimensional character of 
the dynamics of flavor symmetry breaking in this problem. More precisely, the situation is the 
following. Using the identity tanh(x) = 1 — 2exp(—2x)/[l + exp(—2x)] and the relation [15]: 

(1 - zy^) exp (-5L.) = f) L«{x)z n , (24) 

where L"(x) are the generalized Laguerre polynomials, the propagator S>E(k) can be decomposed 
over the Landau level poles |H| : 

S E{k) = - iexp (_i) £(-!)■ ^ffif^,, (25) 



with 

D n (eB,k) = (m-fc 3 7 3 



k 2 

1 - irfirf 2 sign(eB))L n (2- L 



eB 
1 + i7i7 2 sifl , n(eS))L n _ 1 (2 



+ 4(^71 + ^272)^-! (2^), (26) 

where L n = L° n and L" 1 (x) = by definition. Then Eq. (P5|) implies that as m — > 0, the condensate 
appears due to the lowest Landau level: 

We would like to note that in 3 + 1 dimensions, the dynamics at the lowest Landau level, described 
by two continuous variables k^ and k^, is two dimensional^. Actually, the lowest Landau level pole 
in S>E{k) is in 3 + 1 dimensions: 

s {k) = ~ texp rw\ k} + B +m * o--™)- (28) 



■"■The fact that a magnetic field reduces the effective dimension of the dynamics of the fermion pairing by two 
units was observed ealier in the theory of superconductivity [M . 



The contribution of this level to the condensate is now: 

0\ \ m^O 



m\0)~-\eB\^(hi^ + O(m u )) "■'-"(.). (29) 



i.e., there is no spontaneous flavor symmetry breaking in a magnetic field in 3 + 1 dimentions. 

In the next section, we will discuss aspects of spontaneous flavor symmetry breaking for (2 + 1)- 
dimentional fermions in a magnetic field in more detail. 

3 More About Dynamical Flavor Symmetry Breaking for 
Fermions in a Magnetic Field 

As was shown in the preceding section, the flavor condensate (0|\I/\I/|0) is nonzero as the fermion 
mass m goes to zero. Although usually this fact is considered as a firm signature of spontaneous 
flavor (or chiral) symmetry breaking, the following questions may arise in this case: 

a) Unlike the conventional spontaneous flavor (chiral) symmetry breaking, the dynamical mass 
of fermions equals zero in this problem. Is spontaneous flavor symmetry breaking "real" in this 
case? 

b) The vacuum |0) was defined as lim m ^ |0) m of the vacuum |0) m in the theory with m ^ 0. 
The vacuum |0) corresponds to a particular filling of the lowest Landau level. Indeed, at m ^ 0, in 
the vacuum |0) m , the states with Eq = m > are empty and the states with Eq = —m are filled, 
i.e. the vacuum |0) = lim m ^ |0) m is annihilated by all the operators a 0p , d 0p and a np , b np , c np , 
d n p (n > 1). On the other hand, at m = 0, there is the infinite degeneracy of the vacuum in this 
problem, connected with different fillings of the lowest Landau level. Why should one choose the 
filling leading to the vacuum |0)? And is there a filling of the lowest Landau level leading to the 
ground state which is invariant under the flavor U(2)7 One might think that the latter possibility 
would imply that spontaneous flavor symmetry breaking can be avoided. 

In this section we will show that there is a genuine realization of the spontaneous breakdown of 
the flavor symmetry in the present problem. More precisely, we shall show that this phenomenon 
satisfies all the criterions of the spontaneous symmetry breaking phenomenon established by Haag 
long ago pifl . We will also discuss such related questions as the status of Nambu-Goldstone (NG) 



modes and induced quantum numbers [22], |23[ in this problem. 



Let us begin by constructing the charge operators Qi = l/2/<i 2 x[\l/^(x),Tj\E'(x)] of the flavor 
U(2) group. By using Eqs.(||), (0) and (|9]), we find: 

oo 

Wo = 2.^1 \ a 0p a ®P ~ "O-p^O-pJ + 2.^1 2-~i \ a np a np ~ np np + C np C np — Ct np Ct np J 
p n=l p 

oo 
W 1 = 1 2-~i I ®0p 0~p <^0— p a 0p ) + ? 2 < 2 < I a np C np C np a np + np u np CL np O np J , 

P 71=1 P 

OO 

V2 = _^ I a 0p^0-p ~r ^O-p^-Op ) T" ^^ ^_, ( a np C np + C np a np + O np a np + d np O np J , I'jUj 

p n=l p 

I „ pi oo 

V3 — ~7j <-5 + _^ ^ a 0p a 0p + "O-p^O-pJ + ^_^ _^ ^ a np a np — ®np®np ~ c n p C np ' "■np^npj 
Z7F p n=l p 

where a np , c np , (6 np , a! np ) are annihilation operators of fermions (antifermions) from the n-th Landau 
level and S = L\Li is the 2-dimensional volume. Now we can construct a set of the degenerate 
vacua: 

|0i,0 2 > = exp(iQ 1 1 +iQ 2 2 )|O) (31) 

where, we recall, the vacuum |0) = lim m ^ |0) m is annihilated by all the operators a np , b np , c np and 
d np . As one can see from Eq. (|30|) , the crucial point for the existence of the continuum set of the 
degenerate vacua is the first sum, over the states at the lowest Landau level, in the charges Qi and 

Q 2 - 

The presence of such a set of the degenerate vacua is a signal of the spontaneous breakdown, 
U(2) — *■ U(l) x U(l). Note that the vacua |0i,02) can be also constructed by replacing the mass 
term m^/^/ by m® Oyfijte eifin where ^fe lt e 2 = exp(zQi0i + z<3202)^, and then performing the limit 
m — > 0. Again, this is a standard way of constructing degenerate vacua in the case of spontaneous 
breakdown of a symmetry. 

One can check that different vacua |0i, 2 ) become orthogonal as size L x in the Xi-direction goes 
to infinity. For example: 

|(O,0 2 |O,0 2 )| =[]|cos0| =exp(L! fdkln\cos9\), = 2 - 2 , 
p J 

and at 7^ or 7r , it goes to zero as L\ — > oo ( the vacuum |O,0 2 + n) = — |O,0 2 )). It also goes 

to zero as the maximum momentum |/c max | = A ( A is an ultraviolet cutoff) goes to infinity. As 
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usual, this point reflects the fact that spontaneous symmetry breaking occurs only in a system with 
an infinite number of degrees of freedom. One can check that in this case all states (and not just 
vacua) from different Fock spaces F{ 8l g 2 }, defined by different vacua \9i,9 2 ), are orthogonal. That 
is, different vacua |#i, #2) define nonequivalent representations of canonical commutation relations. 
On the other hand, taking the ground state 

\n) = cjd f i(9 1 ,9 2 ,9 3 )\9 l ,9 2 ), (32) 

where d/u, is the invariant measure of SU(2) and C is a normalization constant, we are led to the 
vacuum \Q) which is a singlet with respect to the flavor U(2). In fact, the set of the vacua {|#i, #2)} 
can be decomposed in irreducible representations of 577(2): 

{|*i,W = {l^ } >} (33) 

Why should we consider the vacua \9i,9 2 ) instead of the vacua |fiW)? 

To answer to this question, we consider, following Haag []2l[| , the clusterization property of 
Green's functions. It means the following. Let us consider a Green's function 

G (n +k) = (0 | T JJ A(a ..) f[ Bj ( yi )\0), (34) 

i=\ j=l 

where Ai(xi), Bj(jjj) are some local operators. The clusterization property implies that when 
rfj — ► 00 (rfj = (xj — yj) 2 ) for all i and j, the Green's function then factorizes as follows: 

n k 

G (n +k ) _> (0 |Tn^(^)|0)(0|Tn^(%)|0). (35) 

i=l 3=1 

The physical meaning of this property is clear: clusterization implies the absence of instantaneous 
long-range correlations in the system, so that the dynamics in two distant spatially-separated 
regions are independent. 

The clusterization property takes place for all the vacua | (9i , 6*2) - The simplest way to show this 
is to note that the vacuum \9i,9 2 ) appears in the limit m — > from the vacuum in the system with 
the mass term m^> 6x62^ &i8i- Since at m 7^ 0, the vacuum in this system is unique, the clusterization 
is valid at every value of m 7^ 0. Therefore it is also valid in the limit m — > 0, as far as the Green's 
functions exist in this limit. In connection with that, we would like to note that, in thermodynamic 
limit Li,L 2 — > 00, the vacuum ^1,^2) is the only normalizable and translation invariant state 

11 



in the Fock space Fq x q 2 . To show this, let us introduce the operators a n (k) = (Li/2ir) 1 ' 2 a 



np- 



b n (k) = (L^tt) 1 ^ c n (k) = (L 1 /2tt) 1 / 2 c„ p , d n {k) = (L.^nY^dnp, where k = 2<k P /L x . They 
satisfy the commutation relations [a n (k), a n ,(k')} = 5 nn >5(k — k'), etc. Therefore, though states of the 
form TTj %{kj) T7j do{kj)\9i, #2) have zero energy, they are not normalizable and, at J2i h + J2j kj ¥" 0> 
not translation invariant. 

On the other hand, the clusterization property is valid not for all Green's functions on the vacua 
|fiw). As an example, consider the Green's function 

G (4) = (fi|T(*(xi)*(x 2 ))(*(yi)*(y 2 ))|fi), (36) 

where \Q) is the vacuum singlet (|32l). Since the bilocal operator ^/(xi)^/(x 2 ) is assigned to the 
triplet of SU(2), the clusterization property would imply that 

G W _> (n\T^(x 1 Mx 2 ))\Q)(Q\m yi My 2 ))\Q) - 0. (37) 

as rfj = (x, ; — y^) 2 — > 00. However, since 

(ft|T($(a;i)tf(x 2 ))|ft< 3 >) ^ 0, 

(OWlT^CyO*^))^) ^ 0, (38) 

where |f2^ 3 ' ) ) is a state from the vacuum triplet, we see that G^ does not vanish as rf- — ► 00. 

Thus the clusterization property does not take place for the |f2^)-vacua. 

This is a common feature of the systems with spontaneous continuous symmetry breaking j|, pTf : 
an orthogonal set of vacua can either be labeled by the continuous parameters {#«}, connected with 
the generators Qi of the broken symmetry, or it can be decomposed in irreducible representations 
of the initial group. However, the latter vacua do not satisfy the clusterization property. 

All the Fock spaces F^g^y yield physically equivalent descriptions of the dynamics: in the 
space -F{0 1( 9 2 }, the SU(2) spontaneously breaks down to U{e 1 e 2 }iX)i where the ^7{e 1 e 2 }(l) symmetry 
is connected with the generator Q\ 1 2 = exp(iQi6i + iQ 2 8 2 )Q 3 exp(—iQi6i — iQ 2 9 2 ). Are there 
NG modes in the present system? To answer to this question, let us consider the thermodynamic 
limit Li,L 2 —* 00. One can see that in every Fock space F^^y, with the vacuum \9i,9 2 ), there 
are a lot of "excitations" with nonzero momentum k and zero energy E created by the operators 
a Q (k) and d (k). However there are no genuine ( i.e. with a nontrivial dispersion law) NG modes: 
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the energy E is E = at the lowest Landau level. Since the Lorentz symmetry is broken by a 
magnetic field, this point does not contradict to Goldstone's theorem^. This of course does not 
imply that the existence of NG modes is incompatible with a magnetic field: the situation is model 
dependent. As will be shown in Sees. 4-7, even the weakest attractive interaction in the problem of 
2 + 1 dimensional fermions in a magnetic field is enough to "resurrect" the genuine NG modes. The 
key point for their existence is that the flavor condensate (0|\E , \E'|0) and the NG modes are neutral, 
and the translation symmetry in neutral channels is not violated by a magnetic field (see the next 
section). We shall also see that the "excitations" from the lowest Landau level (with quantum 
numbers of the NG modes) in the problem of free fermions in a magnetic field can be interpreted as 
"remnants" of the genuine NG modes in the limit when the interaction between fermions is being 
switched off. Moreover, we shall see in Sec. 5 that the vacua {61,62) constructed above yield a very 
good approximation for the vacua of systems with weakly interacting fermions in a magnetic field 
(in fact, it appears that the role of the vacua {61,62) is the same as that of the 9- vacua of the ideal 
Bose gas for an almost ideal Bose gas in the theory of superfluidity ||). 



In conclusion, let us discuss the phenomenon of induced quantum numbers |22|, |23| in this 
problem. As it follows from Eg. flSTf) , the vacuum {61,62) is an eigenstate of the density operator 
p\ 1 = Mms^ocQi 1 /S with a nonzero value: 

pl ei62} \6 1 ,6 2 )= l -^\6 1 ,6 2 ). (39) 

Thus, there is the induced quantum number of the operator p\ 1 in the \6\62) vacuum^. This 
fact is intimately connected with the phenomenon of spontaneous flavor symmetry breaking in 
this problem. Indeed, since Q\ x is one of the generators of the non-Abelian SU(2) symmetry, 
its vacuum eigenvalue would be equal to zero if the symmetry were exact and the vacuum were 
assigned to the singlet (trivial) representation of SU(2). This is in contrast to the case of Abelian 
£7(1) symmetry: since U(l) has an infinite number of one dimensional representations, the vacuum 
can be an eigenstate of the charge density p = lim^oo Q/S with an arbitrary eigenvalue in that 
case. 



2 As to a nonrelativistic analog of the Goldstone theorem, it has been proved only for translation invariant systems 
with finite range interactions |24fl . Systems in a magnetic field do not satisfy this condition. 
3 Note that this fact agrees with the consideration in Ref. 123]. 
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Note that, since the SU(2) is spontaneously broken here, it is appropriate to redefine the gen- 
erator of the exact U^g^l) symmetry as Ql 1 = Q\ 1 — \eB\S/2ix. 



4 The Nambu-Jona-Lasinio Model in a Magnetic Field. 
General Consideration. 

In this and the following four sections, we shall consider the NJL model in 2 + 1 dimensions. This 
model gives a clear illustration of the general fact that a constant magnetic field is a strong catalyst 
of generating a fermion dynamical mass in 2 + 1 dimensions. 

Let us consider the 2 + 1 dimensional NJL model invariant under the U(2) flavor transformations: 



£ 



tf,Z7"A»* 



G 

+ 2 



(M) 2 + (^ 7 5#) 2 + 



+ (*73*) 5 



(40) 



where D^ is the covariant derivative (0) and fermion fields carry an additional, "color", index 
a = 1,2, . . . ,N. This theory is equivalent to a theory with the Lagrangian density 



C 



ll^i-fD^ 



\& (a + 7 3 r + ry 5 7r) \& 



2G 



a + ix + r 



(41) 



The Euler-Lagrange equations for the auxiliary fields a, r and tx take the form of constraints: 



-G(W),t= -G(^7 3 ^),vr = -G(^7 5 ^). 



(42) 



The Lagrangian density fl4i"D reproduces Eq. (140|) upon application of the constraints ( 42]) . 

The effective action for the composite fields is expressed through the path integral over fermions: 

1 



T(o,t,ix) 



exp (if ) 



2G 



d 6 x{o l + -r + 7T Z ) + Via, T, TX 



5 ' J ' '/) 



[d\I>] [d\P] exp <{ - / <i 3 x 



(a + 7 3 r + Z7 5 7r) I* 



(a + 7 3 r + Z7 5 7r) 



if/Ai^D. 



exp < TrLn 



WDn - 



(43) 



(44) 
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i.e.. 

„3„ , • ,5 



r(er, r, n) = -iTrLn i^D^ - (a + Y T + ry 5 7r) . (45) 

As N — > oo, the path integral over the composite (auxiliary) fields is dominated by stationary 
points of the action: 5T/5a = 5T/5t = 5T/5ti = 0. We will analyze the dynamics in this limit by 
using the expansion of the action T in powers of derivatives of the composite fields. 

Is the 1/N expansion reliable in this problem? This question appears naturally since, as was 
emphasized in Sec. 2, a magnetic field reduces the dimension of the dynamics of the fermion pairing 
by two units. If such a reduction took place for the whole dynamics (and not just for that of the 
fermion pairing), the 1/N perturbative expansion would be unreliable. In particular the contribution 
of the NG modes in the gap equation, in next-to-leading order in 1/N, would lead to infrared 
divergences. Just such a situation takes place in the 1 + 1 dimensional Gross-Neveu model with a 
continuous chiral symmetry [BSJ. This phenomenon reflects the Mermin- Wagner- Coleman theorem 



forbidding spontaneous breakdown of continuous symmetries in space dimensions lower than 
two. 

Fortunately, as will be shown in Appendix C, this is not the case in the present problem. The 
central point is that condensate (Ol^^lO) and the NG modes are neutral in this problem. As we 
shall see in Sec. 6 and Appendix C, this is reflected in the structure of the propagator of the NG 
modes : unlike the fermion propagator, it has a genuine 2 + 1 dimensional structure. As a result, 
their contribution to the dynamics does not lead to infrared divergences, and the 1/N expansion is 
reliable in this problem. This point is intimately connected with the status of the space-translation 
symmetry in a constant magnetic field. In the gauge (Q), the translation symmetry along the 
a^-direction is broken (though it can be restored by applying a certain gauge transformation). 
Therefore the momentum hi is a bad quantum number for fermions and all other charged states ( 
see Eqs. @ and (0)). However, for neutral states , both the momenta k\ and k 2 of their center of 
mass are conserved quantum numbers (this property is gauge-invariant) p7fl . In order to show this 
fact in the gauge (H), let us introduce the following operators describing space translations in first 
quantized theory: 

P X1 = \^~ , P X2 = -jT- + QBx u (46) 

l OX 1 l ox 2 

where Q is the charge operator. One can easily check that these operators commute with the 
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Hamiltonian of the Dirac equation in a constant magnetic field. Also, the commutator [P X1 , P X2 ] i s: 

[P xl ,P xa ] = -iQB. (47) 

Therefore the commutator equals zero for neutral states, and both the momenta k\ and k 2 can be 
used to describe the dynamics of the center of mass of neutral states. As we shall see, this point is 
important for providing the (2 + l)-dimensional character of this dynamics. 

5 The NJL Model in a Magnetic Field. The Effective Po- 
tential. 

We begin the calculation of T by calculating the effective potential V. Since V depends only on the 
S'C/(2)-invariant p 2 = a 2 + r 2 + 7r 2 , it is sufficient to consider a configuration with r = it = and a 
independent of x. So now T(a) is 

f(a) = -iTrLn(iD - a) = -iLnDet(iD - a), (48) 

where D = 7 M -D M . Since 

Det(iD -a) = Det(-f 5 (iD - a)rf) = Det(-iD - a), (49) 

we find that 



T(a) = — Tr 
K ' 2 



Ln(W -a) + Ln(-iD - a) 



= --TrLn{D 2 + a 2 ). (50) 

Therefore T(a) can be expressed through the following integral over the proper time s: 

T(a) = --TrLn(D 2 + a 2 ) = - f d 3 x f°° —tr(x\e- ia ^ 2+ ^\x) (51) 

2 2 J Jo s 



where 



D 2 = D^--^YF^= (52) 
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The matrix element (x\e %s<yD +CT )\y) can be calculated by using the Schwinger approach |T7|]. It is: 



[x\e 



where 



i S (D 2 +<T 2 )ly\ 


— 


e- lS(72 (x\e- lsD ^\y) 




+ 


7 X 7 2 sin(eBs) 


= 




= 


-i- 
e 4 .-Use 


2 -Scl) . 




8(tts) 3 / 2 





cos(eBs) + 



eBscot(eBs) +77 eBs 



rx 1 f (f 2 y* 

S d = e\ y Ardz x --(x-y)Jg^+ [ -^- 



[1 — eBs cot(eBs)} (x — y)^. 



(53) 



(54) 



Here the integral /* A^dz is taken along the straight line. 
Substituting Eq. (B^) in Eq. fl5"T|), we find 



IV) = ^L^ /> x r ^e- isa2 eBscot{eBs). 



Ait 3 / 2 J Jo s 

Therefore the effective potential is 

vv) = L + K( P ) = £ + J^ /; ^e-V eBscoth(eBs ), 

where p 2 = a 2 + r 2 + 7r 2 , and now we introduced explicitly the ultraviolet cutoff A. 
By using the integral representation for the generalized Riemann zeta function £ 



/ cis^- 1 e- /3s coths = r(/i) 
Jo 



a^Cte §) - /T" 



(55) 



(56) 



(57) 



which is valid at \i > 1, and analytically continuing this representation to \i — —\, we can rewrite 



Eq. (69) as 



vw = N - 

TX 

_P_ 
2/ 2 



+ 0(1/ A), 



A./ 1 (p/) : 



2' 2 



+ 1 



(5f 



where the magnetic length / is I = \eB\ 1//2 and here we introduced the dimensionless coupling 
constant 

g = N-G. (59) 

7T 
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(We recall that C(/i, |) is defined as 



<^t^W m 



at /i> 1 |T§.) 

Let us now analyze the gap equation dV/dp = 0. It is 

Ap p f°° ds 



ng 27r 3 / 2 Ji/A 2 s 3 / 2 
which can be rewritten as 

„2;2 



/ -^e~ sp eBscoth(eBs), (61) 

Jl/A 2 S 6 '^ 



2A/(i-4=)P = y + v / 2pcdl + ^-) + 0(l/A). (62) 



As S — > 0, we recover the known gap equation |14 



P 2 = M ( 4 - - ) ■ (63) 



^ 9, 

It admits a nontrivial solution only if the coupling constant g is supercritical, g > g c = ^pn (as 
Eq. (|4T|) implies, a solution to the gap equation, a = a, coincides with the fermion dynamical mass, 
& — m dyn)- We will show that the magnetic field changes the situation dramatically: at B ^ 0, a 
nontrivial solution exists at all g > 0. The reason for this is that the magnetic field enhances the 
interaction in the infrared region (large s): at B ^ 0, the integral in Eq. ([Jl]) becomes proportional 
to 1/p as p — > 0. 

We shall first consider the case of subcritical g, g < g c = ^/k, which in turn can be divided into 
two subcases: a) g <C g c and b) g ^ g c — (nearcritical g). Assuming that \al\ <C 1 at g <C g c , we 

find from Eq. fl62"|): 

\eB\g^Tx 

-*■ s " K 2Afe^y- (64) 

Since Eq. (|6lD implies that the condition |aZ| <C 1 fulfills at all (7 satisfying (g c — g) ^$> ^-4 — , the 
relation (|^) is actually valid in that whole region. 

Note the following interesting point. Eq.(g2|) implies that m dyn = (0|er|0) = — 7r^/A r A(0| 1 i r \I/|0>. 
From here and Eq.(H) we find that the condensate (Ol^lO) is (0|W|0) = -N\eB\/2ir in leading 
order in g, i.e. it coincides with the value of the condensate calculated in the problem of free 
fermions in a magnetic field (see Eq. (|22|) ) . This point implies that at small g (weakly interacting 
fermions) the \$i, #2)~vacua constructed in Sec. 3 are good trial states for the vacua of the problem 



with interacting fermions. This point also explains why the dynamical mass md yn i n this problem 
is an analytic function of g at g = 0: indeed, the condensate exists already at g = 0! 



At g c — g<J\eB\/A, introducing the scale m* = A(l/g — l/g c ), we get the equation 

2m ' l -Wi + Vk i^ + 1 ) (65) 



which implies that in the nearcritical region, m dyn is 



m dyn = o ~ \eB\ 1/2 . 



(66) 



Thus in the scaling region, with g c — g<J\eB\/A, the cutoff disappears from the observable m^. 
This agrees with the well-known fact that the critical value g c = ypK is an ultraviolet stable fixed 
point at leading order in 1/N ||14| . The relation ( |66| ) can be considered as a scaling law in the 
scaling region. 

In the supercritical region, at g > g c , the analytic expression for m^,, can be obtained at weak 



\eB\, satisfying the condition J\eB\/m dyn <C 1, where m d y n is the solution of the gap equation ( |6~^ 



with B = 0. Then, using the asymptotic formula |18 

1 



£(z, q)q -»• 06 



(z - Vjq*- 1 



we find from Eq. (|62] 



m d . 



u" 



a = m 



(o) 

dyn 



1 + 



{eBf 



12 (m 



(0) N 



(67) 



(6? 



i.e., m^n increases with S. The numerical study of the equation ( |6^) shows that m dyn increases 
with B at all values of g and B. 

A striking fact is that, unlike the gap equation ( |B"3|) with B = 0, the gap equation with 5^0 
does not have the trivial solution a = 0. Indeed, Eq. (|56[) implies that <i\^/cicr| cr=0 = dV/da\ a= o, 
and then we find from Eqs. (|4TD and (0) that 

U Dl 

(69) 



(/K U =0 = (0|M|0)| 9=0 = -A^l ^ 



(icr 



2^ 



(see Eq. (|22|) ). Thus, despite the spontaneous character of the U(2) symmetry breakdown, there is 
no trivial solution (stable or unstable) in the magnetic field at all values of g. 
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6 The NJL Model in a Magnetic Field. The Kinetic Term 
in the Effective Action. 



Let us now consider the kinetic term Ck in the effective action fl43|). 
The U{2) symmetry implies that the general form of Ck is 



A = N^-{d,p 3 d uPj ) + N^-{ P] d,p 3 ) ■ 
z p 

(p l d, Pt ) (70) 



where p = (cr, r, n) and F^" , F% v are functions of p 2 = a 2 + t 2 + 7T 2 . To find the functions F^ v , F% v , 



one can use different methods. We used the method of Ref . p8[ . The derivation of Ck is considered 
in Appendix A. Here we shall present the final results. 

The functions Ff and Ff take the form F^ = g^Fj*, Ff = g^Ff where 



§tt VV2 V 2' 2 

F ll = F 22 



F 2 °° 



47rp 
Z 



(^4^-M<H 



t2 ~ t2 ~8n 



16vr \2\/2 v 2' 2 



V2 



+ v / 2(p/) 2 C(^^ + l)+2p/-(p/)- 1 

(we recall that the magnetic length / = \eB\~ l l 2 ). 

We would like to emphasize that, as follows from Eq.(|7T|), the propagator of the NG modes in 
leading order in 1/N has a genuine (2 + l)-dimensional form. We shall see in Appendix C that this 
fact is crucial for providing the reliability of the 1/N expansion in this problem (physical reasons 
for the (2 + l)-dimensional character of the dynamics of the neutral NG bosons are considered in 
the next section). 

Now, knowing the effective potential and the kinetic term, we can define the energy spectrum 
(dispersion law) of the collective excitations o and r, it. 
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7 The NJL Model in a Magnetic Field. The Spectrum of 
the Collective Excitations. 

We begin by considering the spectrum of the collective excitations in the subcritical, g < g c , region. 



At g c ~ g > y/\eB\/A (where \al\ < 1, see Eq. (||)), we find from Eqs. (|70|) and fl7T|) the 
dispersion law for the r and n NG (gapless) modes: 

E T , c V2(al)(kY /2 = ^^) (fe 2 ) 1/2 (72) 



(see Eq. fl64|)). As the interaction is switched off, g ^ 0, their velocity, v = gg c (g c — g)~ 1 J\eB\/2A 2 , 
becomes zero, and we return to the dynamics with spontaneous flavor symmetry breaking but 
without genuine NG modes discussed in Section 3. 

In order to define the "mass" (energy gap) M a of the cr-mode, we note that 
d 2 V, „a 2 l r°° 



da 2 ,a=a 

2tt 



N—r-pr \ dsy/sexp(— (a l) 2 s) coths 

7l 6 ' 2 JO 

' 1 . ,3 (a/) 2 n , A 



V2 v 2' 2 
(see Eq. (||)). Then we find from Eqs. fl70|), (J7T|) and (]7§ that 



(73) 



^ 8v^(g c -g) A|eg| i/ 2 (74) 

99cC{\) 

at g c — g ^> J\eB\/A. As g — > 0, the cr-mode decouples (M ff —> oo). 

Thus the dynamics in the problem of a relativistic fermion in an external magnetic field emerges 
from this model in the limit when the interaction between fermions is switched off. The attractive 
(g > 0) interaction "resurrects" the NG modes and they acquire a velocity v ~ g. 



Let us now consider the nearcritical region with g c — g<J\eB\/A. From Eqs. (f70|), (|7I|), we find 
that 

E^ = f{al){k 2 ) 1 ' 2 , (75) 

where 

^=^H^ {j r-^ +{9 ^T ■ (76) 

Since in this nearcritical (scaling) region the parameter a is a ~ |e-B | 1 ^ 2 = l~ x , we conclude that 
the cutoff A disappears from the observables E T and E^ in the scaling region. 
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In the same way, we find from Eqs. (|70|), (|7T| ) and ( |T3[ ) that 

^ ~ \ eB \ (77) 

in the scaling region. 

Let us turn to the supercritical region with g > g c . The analytic expressions for E T ^ and M% 
can be obtained for small \eB\ = l~ 2 , satisfying the condition \al\ ^> 1. Then, using the asymptotic 
formula (^) for zeta functions, we find from Eqs. (0), ([H]) and (|73|): 

E - = t 1 " 8(M (fe2,1/2 ' (78) 

M * - « 2 (> - w) • < 79 > 

where <r is given in Eq. ( j6Jj|) . These relations show that the magnetic field leads to decreasing both 
the velocity of the NG modes (it becomes less than 1) and the mass (energy gap) of the a mode. 

Let us indicate the following interesting point intimately connected with the (2 + l)-dimensional 
character of the dynamics of the neutral NG modes. The (2 + l)-dimensional character is reflecded 
in that the velocity v„- )T = dE n ^ T /dk is not zero. As follows from Eqs.(|72|), (ff5|) and (|78D , the velocity 
v n>T decreases with va^ yn — & and becomes zero (i.e. the dynamics becomes (0 + l)-dimensional) 
when rridyn — > 0, i.e. when the interaction is switched off (g -^ 0). The reason for this is clear: 
since at g — the energy of the neutral system made up of a fermion and an antifermion from 
the lowest Landau level is identically zero, its velocity is also zero. This fact in turn reflects the 
point that the motion of charged fermions in the a;i-X2 plane is restricted by a magnetic field. On 
the other hand, at g > 0, there are genuine neutral NG bound states (with the bound energy 
AE WtT = 2md yn — E ntT \k=o = Zmdyn)- Since the motion of the center of the mass of neutral bound 
states is not restricted by a magnetic field, their dynamics is (2 + l)-dimensional. 

Let us now discuss the continuum limit A — > oo in more detail. As is known, at B — 0, in 
this model, an interacting continuum theory appears only at the critical value g = g c = -*/¥ (the 
continuum theory is trivial at g < g c ) ||, [TJ], [15], 16| . Therefore, since at g < g c , in the continuum 



limit, there is no attractive interaction between fermions, it is not surprising that at g < g c , the 
dynamical mass m^n ~ g\eB\/A disappears as A — >• oo. 

At B = 0, the continuum theory is in the symmetric phase at g ^ g c — and in the broken 
phase at g — > g c + . On the other hand, as follows from our analysis, in a magnetic field , it is in 
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the broken phase both at g — > g c — and g — ► g c + (though the dispersion relations for fermions 
and collective excitations p are different at g ^ g c — and g — > g c + ). 

Up to now we have considered four-component fermions. In the case of two-component fermions, 
the effective potential, V 2 is V 2 (cr) = V{o)/2 where V{o) is defined in Eqs. fl56|) and fl58l). However, 
the essential new point is that there is no continuous (U(2)) symmetry (and therefore NG modes) in 
this case. As in the case of four-component fermions, in an external magnetic field, the dynamical 
fermion mass (now breaking parity) is generated at any positive value of the coupling constant g. 

The NJL model illustrates the general phenomenon in 2 + 1 dimensions: in the infrared region, 
a magnetic field reduces the dynamics of fermion pairing to a one-dimensional dynamics (at the 
lowest Landau level), thus catalysing the generation of a dynamical mass for fermions. A concrete 
sample of dynamical symmetry breaking is of course different in different models. 



8 Thermodynamic Properties of the NJL Model in a Mag- 
netic Field. 

In this section, we will study the thermodynamic properties of the NJL model in a magnetic field. 
In particular, we will show that there is a symmetry restoring phase transition at high temperature. 
Our goal is to determine the thermodynamic (effective) potential in the NJL model in a magnetic 
field. Although we are mostly interested in studying the system at finite temperature T and zero 
chemical potential \x (i.e., at equal densities of fermions and antifermions), we shall derive the 
effective potential Vp^(a) (at the leading order in 1/iV) at arbitrary values of /3 = 1/T and \x. 

Since in the leading order in 1/N, the effective potential V(a) = Vp ttl (a)\g _ is given by a 

/i = 
sum of one- (fermion) loop diagrams, the thermodynamic potential is in this approximation: 



VeM 



V(a) + V^(a) 



N 

71 

- N- 
+ e 



A 



2A 9 
eB 



-Da' 



v/2\/ 1 (alf 



< 



2n{3 



P ^V 2' 2 

oo 

ln(l + e- /3(ff ^ ) )+2^1n(l + 



fc=i 



- /3( v / ^- ) ) + (/i ^_ /i) ' 



a 
2P 



JO) 
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(see Eq. (|58|)). Here the sum is taken over all the fermion and antifermion (with fi — > —fj,) Landau 
levels; the factor \eB\/2ir describes the degeneracy of each level. The relation fl30|) is derived (in 
the framework of the imaginary time formalism |2J|) in Appendix B. We also show there that it 
can be rewritten as 

- 2 N r™ dt 



VpA°) 



a 



+ 



o 



2G 4vr 3 / 2 / 3 

\2^ Hl 4:7ivr , 



t 3/2 



e -^ 2 )coth£- 



$1} 



where 64 is the fourth Jacobian theta function |fLq| . 

Henceforth we will consider the case of zero chemical potential corresponding to equal densities 
of fermions and antifermions in the system. In this case, the thermodynamic potential Vp = V^ iAt | M= o 

is: 

A 



V {ar) 



N 

71 

a 
2P 



2^ 9 



Da' 



7T{3 
00 

+ 2]Tln(i + e~ 



ln(l + 6-1*) + 



/ 3 ^ 



1 {all 
2' 2 



2 1 2k 



Ps ° 2 +% 



fc=l 



(82) 



We solved numerically the gap equation 



dVpja 

da 



0. 



J3) 



The main result is that at T = T c ~ md yn (actually at T = T c ~ Tridyn/2), there is a symmetry 
restoring (second order) phase transition (see Figs. 1 and 2). The phase diagram in the B — T plane 
is shown in Figs. 3 and 4. 

We recall that there cannot be spontaneous breakdown of a continuous symmetry at finite 
(T > 0) temperature in 2 + 1 dimensions (the Mermin- Wagner- Coleman (MWC) theorem ||26|| ) . 
This happens because at nonzero temperature the dynamics of a zero mode in 2 + 1 dimensional 
field theories is two dimensional. As a result, strong fluctuations of would be NG modes lead to 
vanishing the order parameter connected with a spontaneous breakdown of a continuous symmetry. 
In the NJL model with a finite temperature (both at B = and in a magnetic field), the MWC 
theorem manifests itself only beyond the leading order in 1/N. One plausible possibility of what 
happens at T 7^ beyond the leading order in 1/N is the following. The dynamics of the zero mode 
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in this model is essentially equivalent to that of the SU(2) cr-model in 2 dimensional Euclidean 
space. As is known, the SU(2) symmetry is exact in the latter model and, as a result, the would 



be NG bosons become massive exitations |30 ]. Therefore it seems plausible that in the (2 + 1)- 
dimensional NJL model in a magnetic field, the SU(2) symmetry will be restored at any finite 
temperature, and the dynamically generated mass m^ n of fermions will disappear. 

The question whether this, or another, scenario is realized at finite temperature in this model 
deserves further study. 

9 Conclusion 

The main result of this paper is that a magnetic field is a strong catalyst, generating a fermion 
mass (energy gap), in 2 + 1 dimensions. It would be worth considering the present effect in 2 + 1 
dimensional effective theories describing high temperature superconductivity and the quantum Hall 
effect where a magnetic field is an important ingredient of the dynamics. In connection with this, 
we note that in some models of high temperature superconductivity of Ref. 0, the energy gap 
in the electron spectrum results from electron-hole (i.e., fermion-antifermion rather than fermion- 
fermion) pairing. Also, using the four-component spinors in these models reflects the presence of 
two sublattices in high temperature superconductors. Another, potentially interesting, application 
of the present effect may be in 3 + 1 dimensional field theories at high temperature. Since at 
high temperature, their dynamics effectively reduces to that of 2 + 1 dimensional theories, it might 
happen that in a magnetic field, at high temperature, fermions (quarks in quark-gluon plasma, for 
example) acquire a dynamical mass and NG excitations appear. 

It would be interesting to check the realization of this effect in 2 + 1 dimensional theories in 
lattice computer simulations. Note that the recent computer simulations of the 2 + 1 dimensional 



NJL model |L6| show that the 1/N expansion is quite reliable, at least at iV > 12. 

The essence of the present effect is that in a constant magnetic field, the dynamics of fermion 
pairing is one-dimensional: the pairing takes place essentially for fermions at the (degenerate) lowest 
Landau level. This implies the universal character of this effect in 2 + 1 dimensions. 

In this paper, we considered the dynamics in the presence of a constant magnetic field only. It 
would be interesting to extend this analysis to the case of inhomogeneous electromagnetic fields in 
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2 + 1 dimensions. In connection with this we note that the present effect is intimately connected 
with the fact that in 2 + 1 dimensions, the massless Dirac equation in a constant magnetic field 
admits an infinite number of normalized solutions with E — (zero modes); more precisely, the 
density of such solutions is finite. One may expect that the same effect will take place for any 
electromagnetic field configuration in which the density of zero modes is finite. As we have known 
recently, the program of the derivation of a low energy effective action in 2 + 1 dimensional QED 
in external electromagnetic fields has been developed in Ref. [|3l] . 
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Appendix A 

In this Appendix, we derive the expressions for the fermion propagator and for the kinetic term 
Ck in the effective action. 

In the coordinate space, the fermion propagator is 

S{x,y) = {iD + m) x (x\ — ^-\y) = 

= (iD + m) x ds(x\exp\-is(m 2 + D 2 )]\y), (84) 

Jo 

where D = 7^-0^ and D^ is the covariant derivative in Eq. ([|). 

The matrix element (x\e~ ls ^ m +D ^\y) can be calculated by using the Schwinger (proper time) 

approach [O. It is 

— t- 

(x\e- ls{m2+f)2) \y) = * * / ^ - sm2 HeBscot{eBs)+ (85) 

Sins) 6 ' 2 v 

+ 7 1 7 2 e5s 
where 

Sd = ej A e fdz x - -(x - y) v [<f» + U J } (l 
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- eBs cot(eBs)) (x - j/) M . (86) 

Here the integral is calculated along the straight line. 

From Eqs. fl84] ) and (|85|) , we find the expression (|T8| ) for the fermion propagator. 

Let us now consider the derivation of the kinetic term (70) in the low-energy effective action: 



A = N-^-{d l ,p j d v p j ) + N-^-{pjd ll pj) 



Ud, Pj d vPj ) + n^ 2 

z p 

{pidvpi), (87) 



where p = (a, r, n) and F^ , F% v depend on the C/(2)-invariant p 2 = o 2 + r 2 + it 2 . The definition 
Tfc = J d 3 xCk and Eq. (|87| ) imply that the form of the functions F^", F^ is determined from the 
equations: 



ST 

N JZ7ZAxZ7a\\cr = const = ~( F i + 2F 2 )la = const ' d ^$ 



5a(x)5a(0) 



x 



T = IT = T = 7T = 



# 2 r 

A-7r('^U-7r('nV cr = const = ~ F ± I a = const ' "»°vO \x). (89) 



Here T^ is the part of the effective action (|43|) containing terms with two derivatives. Eq. ( fJ3] 
implies that T^ = r\. Therefore we find from Eq. (|89|) that 



Ff = -— [d 3 xx»x % /' r * . = -— fd 3 xxV- 
1 2 7 57r(x)(57r(0) 2 J 

<5 2 f 

(90) 



57r(x)57r(0) 

(henceforth we shall not write explicitly the condition r = it = 0, cr = const). Taking into account 
the definition of the fermion propagator, 

iS- 1 =ib-a, (91) 



we find from Eq. (f45|) that 

7 , ? 77^ = -^r(5(x, 0)^7^(0, x)^ 7 E 
= -itr(5(x,0)i 7 5 ^(0,x)27 E 
= -i J < ^^e^tr{S{k)i 1 "~S{k + q^) (92) 

(the functions S(:r) and S(k) are given in Eqs. (p^8|)-(pOD). Therefore, 

itf" 1 /■ d 3 k (~ 5 d 2 S(k) . A 

Fr =-—J w? tr ( 5( ^ 7 ^c 7 J • (93) 
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In the same way, we find that 

iN- 1 



Ff 



4 

iN- 1 



** tr(S(k)?M) ~ W 



(2tt) 3 V v 'dkMv J 2 



v^\jnj V 



(94) 



d 3 k ( \,.,d 2 S{k) ~ n _^ b d 2 S{k) . 5 N 



; tr S(k) 



dk^dk v 



- S(k)ir 



dk^dK 



t'y 



4 J (2tt) 3 

Taking into account the expression for S(k) in Eq. (p0|) (with m = a), we get: 
d 2 S(k) 



dk°dk° 



d 2 S(k) 
WdU 



/■OO 

= 2il" / dttexv{(R(t)){a(l + rrfi 2 T) + 
Jo 

+ 3A; 7°(1 + thVT) - fey (1 + T ') + 2«t/ 2 (A; )V(l + 

+ thVT) + 2^/ 2 (A; ) 3 7 (l + ^V^) - 2ttl 2 (k°) 2 (k i 1 i )(l + T 2 )}, 

= -2il 4 / dtTexp( J R(t)){cr(l + r77 1 7 2 T)- 

JO 

- fcY(l + ^ 2 ) - 2ifc J V(l + T 2 ) + fe °7°(l + Vl^T) - 

- 2iTl 2 (W) 2 a(l + rjj^T) - 2iTl 2 (W) 2 k° 1 °(l + 
+ ?77 1 7 2 r) + 2iTl 2 {k j ) 2 k i Y{l + T 2 )} 



(95) 



(96) 



(i,j = 1, 2; there is no summation over j), where 

rj = sign(eB), T = tant, 
R(t) = -tt(al) 2 + it(k°) 2 -il 2 k 2 T (97) 

Eqs. (|3)' dUD' (0) i m ply that nondiagonal terms F^ v and Fjf are equal to zero. The diagonal 



terms are determined from Eqs. (|20|) , (|93| ) -(|96 D , after rather long, although routine, calculations: 
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(101) 
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Here (besides Eq. (|57|)) the following relations were used [p"8| : 
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dr = 2 1 -T(/*)[2C(ji-l,£ 
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Appendix B 



In this Appendix we shall derive the thermodynamic potential Vp^ in the NJL model 
(3 = 1/T is an inverse temperature and /i is a chemical potential. 
As is well known [p9f| , the partition function 



(102) 
(103) 

(104) 
; here 



'Aa< 



jy 



exp (-/?#'; 



(105) 



is expressed through a path integral over fields of a system (here H' = H — it J \I/7 \I/d 2 £, H is the 
Hamiltonian of the system). In the NJL model ([i"0"D, (f4T|), the path integral is: 

-i/3 

/o 



^ 



[dtf] [d#] [da] [dr] [dvr] exp {i / dt d 2 x 
1 



* S "* " 2G" 



} 



(106) 



where p 2 = cr 2 + r 2 + ir 2 , S is the fermion propagator ( |i~8|) with m replaced by a + 7 3 r + ry 5 7r, and 
while the fermion fields satisfy the antiperiodic boundary conditions, 



* 



4=0 



-*!*=_ 



t=-i@i 



* 



i=0 



-tfl 



t=-i(3, 



(107) 



the boson fields satisfy the periodic boundary conditions. 

In order to calculate the thermodynamic potential Vp }ll (p), it is sufficient to consider configura- 
tions with t = 71 = and a = const. Then the potential is defined as: 



-i(3 



"x ■ 



exp{-pV / 3, li [f d 2 x\) = f[d^}[d^}exp{t f % dt j d 2 

}. (108) 

At the leading order in 1/N, this potential defines the thermodynamic properties of the system. 



*«-'* - y> 
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As is known [2?J, in the formalism of the imaginary time, the thermodynamic potential Vp^ can 
be obtained from the representation for the effective potential V, at T = and \x = 0, by replacing 



fc° -> ifj n + //; w„ = ^(2n+l) (109) 



tu n = |(2n + 1) follows from the antiperiodic conditions ( |107| )). Then, using the representation for 



the effective potential in Section 4 and the expression ( p0|) for the fermion propagator, we get: 

a 2 N r°° dt 



V„„(a) = — 



72 



2G 2irf3l 

/jl 2 . A . I 



f°° dt , /2(V 2 ,.2-v 

/ — e" w (<T " M } cotht- 
7o t 



■ Q^nt^nt-), (110) 

where 

oo 

Q 2 (u\t) = 2 ]T e^ T(n+ ^ )2 cos ((2n + l)u) (111) 

n=0 

is the second Jacobian theta function ITH . 



By using the identity |18| 



e 2 Hr) = (-)Vv^e 4 (-|--), (ii2) 



r v r r' 



where 

oo 

e 4 (u|r) = l + 2^(-l)V™ 2r cos(2mz) (113) 

ra=l 

is the fourth Jacobian theta function, one can rewrite the relation flllOp as 



VgJa) = ^ + -^7^/ ^z cotht • 



a 2 N r°° dt 



-t/'ff 3 



^(^Isi'T) ') = (114) 



N r°° dt °° 

V{a) + 2^pJo ^cothtE(-ircosh(^n)- 



• exp-(ta¥ + ^), 

where V(cr) is the effective potential (|oT)|). 

Thus we have derived the representation ( |81~D for the thermodynamic potential. Let us show 
that it is equivalent to the representation (p0|). 
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By using the series 

oo 

coth£ = l + 2^e- 2tm , 

m=l 

the expression for V^p = V^p — V in Eq. ( |114j) can be rewritten as 

N Y( l)n COsh (^ n ) L-flm 



(115) 



V, 



(i,/3 



*W£? 



n 



+ 



(116) 



+ 2 J2 exp(-/3crnWl + 



m=l 



2m 

W 2 



Here we also used the relations [118 



/3 



f°° dxx v ~ 1 exp(-" - jx) = 2(^yl 2 K v (2 x ^ht) 
Jo x 7 



#_$(«) = ^0*) 



'2z' 



.V2,,-* 



(117) 



where K,,(z) is a modified Bessel function. 



Since 



we find that 
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(119) 



m=l 



1 + e 



-2(B<7. 1+ 



^ +2e 



-/3cr./l+ 



(ai) 2 



COS 



h(Ai/9)]}. 



It is now easy to check that the expression for the thermodynamic potential Vp )fl = V+Vp^ coincides 
with that in Eq. (f^). 

Appendix C 

In this Appendix we analyze the next-to-leading order in 1/iV expansion in the (2 + 1)- 
dimensional NJL model at zero temperature. Our main goal is to show that the propagator of 
the neutral NG bosons it and r have a (2 + 1) -dimensional structure in this approximation and 



that (unlike the (1 + 1) -dimensional Gross-Neveu model f25|) the 1/N expansion is reliable in this 
model. 

A review of the 1/N expansion in (2 + l)-dimensional four-fermion interaction models can be 
found in Ref.]l4|]. For our purposes, it is sufficient to know that this perturbative expansion is given 
by Feynman diagrams with the vertices and the propagators of fermions and composite particles 



31 



a, 7r and r calculated in leading order in 1/N. In leading order, the fermion propagator is given in 
Eqs.([T8|)-(pTD. As follows from Eq.(il), the Yukawa coupling of fermions with a, r and 7r is gy = 1 



in this approximation. The inverse propagators of a, r and n are [14, 2£ 



(120) 



Dp\x) = N (—5 3 (x)+itr[S(x,0)T p S(0,x)T p ] 

where p = (cr, r, n) and T a = 1, T T = 7 3 , T n = ry 5 . Here S(x, 0) is the fermion propagator ( |18D 
with the mass vn,^ yn = a defined from the gap equation d62|) . For completeness, we write down the 
explicit expression for the Fourier transform of the propagators of the NG bosons: 



D-\k) = D-\k) 



Aipin 



du 



sinhs 



exp[— s(aiy 



( 1 — exp[/?(s,w)] J (a/) 2 coshs H — : 
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121^ 



where 



i2(a,«) = J |(/A; ) 2 ( 



/ 2 k 2 cosh s — cosh su 



■u 



(122) 



2 sinh s 

Actually, for our purposes, we need to know the form of these propagators at small momenta only. 
We find from Eqs.(^), ©: 

4ira 



D T {k)=D^k) 



N 



-f\al)[ki-f\al)^} 



21-1 



(123) 



where 



^-(r^&T 



-1/2 



(a0" 



(124) 



(see Eq.(|7|)). 

The crucial point for us is that, because of the dynamical mass radyni the fermion propagator 
is soft in the infrared region (see Eq . (|25|) ) and that the propagators of the r and 7r (123) have a 
(2 + l)-dimensional form in the infrared region (as follows from Eqs. (|70l) , ( |7lf ) the propagator of a 
has of course also a (2 + l)-dimensional form). 

Let us begin by considering the next-to-leading corrections in the effective potential. The 
diagram which contributes to the effective potential in this order is shown in Fig. 5a. Because 
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of the structure of the propagators pointed out above, there are no infrared divegences in this 
contribution to the potential. (Note that this is in contrast to the Gross-Nevue model: because of a 
(1 + l)-dimensional form of the propagators of the NG bosons, this contribution is logarithmically 
divergent in the infrared region in that model, i.e. the 1/N expansion is unreliable in that case). 
Therefore the diagram in Fig. 5a leads to a finite, 0(1), correction to the potential V (we recall 
that the leading contribution in V is of order N). As a result, at sufficiently large values of N, the 
gap equation in next-to-leading order in 1/N in this model admits a nontrivial solution p ^ 0. 
Since the potential depends only on the radial variable p, the angular variables 9 and ip (p = 
(pcos#, psin0cos<p, psin#sin<p)), connected with the r and n, appear in the effective lagrangian 
only through their derivatives. This in turn implies that the r and tc retain to be gapless NG modes 
in the next-to-leading order in 1/N. 

Let us now consider the next-to-leading corrections to the propagators of these NG modes. First 
of all, note that in a constant magnetic field, the propagator of a neutral local field p(x), D^x, y), 
is translation invariant, i.e. it depends on (x — y). This immediately follows from the fact that the 
operators of space translations (|4q) take the canonical form for neutral fields (the operator of time 
translations is id/dt for both neutral and charged fields in a constant magnetic field). The diagrams 
contributing to the propagators of the NG modes in this order are shown in Fig.5b. Because of the 
dynamical mass m^ n in the fermion propagator, this contribution is analytic at k^ = 0. Since at 
large N the gap equation has a nontrivial solution in this approximation, there is no contribution 
of 0(k°) ~ const in the inverse propagators of r and n. Therefore the first term in the monentum 
expansion of this contribution has the form C\k\ — C^k 2 , where C\ and C<i are functions of al, i.e. 
the propagators take the following form in this approximation: 



D T (k)=D^k) k ="-^f(al) 



1 - i&Ort)) kl - (>(a/) - i<7 2 (<rt)) k 2 ]" 1 . (125) 



(see Eq.fll23D). 

Because of the same reasons, there are also no infrared divergences either in the fermion prop- 
agator (see Fig. 5c) or in the Yukawa vertices (see Fig.5d) in this order. Therefore at sufficiently 
large values of N, the results retain essentially the same as in leading order in 1/N. 

We believe that there should not be principal obstacles to extend this analysis for all orders in 
1/N. 

33 



References 

[1] R. Jackiw, Phys. Rev. D29 (1984) 2375; I. Affleck, Nucl. Phys. B265 (1986) 409; I.V. Krive 
and A.S. Rozhavskii, Usp. Fiz. Nauk 152 (1987) 33 [Sov. Phys. Usp. 30 (1987) 370]. 

[2] A. Kovner and B. Rosenstein, Phys. Rev. B42 (1990) 4748; G.W. Semenoff and L.C.R. Wije- 
wardhana, Phys. Rev. D45 (1992) 1342; N. Dorey and N.E. Mavromatos, Nucl. Phys. B368 
(1992) 614; R. MacKenzie, RK. Panigrahi and S. Sakhi, Phys. Rev. B48 (1993) 3892. 

[3] J. Bardeen, L.N. Cooper and J.R. Schrieffer Phys. Rev. 108 (1957) 1175. 

[4] V.A. Miransky, Dynamical Symmetry Breaking in Quantum Field Theories (World Scientific, 
Singapore, 1993). 

[5] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122 (1961) 345. 

[6] I.V. Krive and S.A. Naftulin, Phys. Rev D46 (1992) 2737. 

[7] A.I. Akheizer and V.B. Berestetsky, Quantum Electrodynamics (Interscience, New York, 1965). 

[8] T. Appelquist, M. Bowick, D. Karabali, and L.C.R. Wijewardhana, Phys. Rev. D33 (1986) 
3704. 

[9] T. Appelquist, D. Nash and L.C.R. Wijewardhana, Phys. Rev. Lett. 60 (1988) 2575. 

[10] E. Dagotto, J.B. Kogut, and A. Kocic, Phys. Rev. Lett. 62 (1989) 1083; Nucl. Phys. B334 

(1990) 279; S. Hands, and J.B. Kogut, Nucl. Phys. B335 (1990) 455. 

[11] D. Atkinson, RW. Johnson, and P. Maris, Phys. Rev. D42 (1990) 602; K.-I. Kondo and H. 
Nakatani, Mod. Phys. Lett. A5 (1990) 407; Prog. Theor. Phys. 87 (1992) 193. 

[12] M.R. Pennington and D. Walsh, Phys. Lett. B253 (1991) 246; R. Pisarski, Phys. Rev. D44 

(1991) 1866. 

[13] M.C. Diamantini, G.W. Semenoff, and P. Sodano, Phys. Rev. Lett. 70 (1993) 3438. 
[14] B. Rosenstein, B.J. Warr and S.H Park, Phys. Rep. 205 (1991) 59. 



34 



[15] K.I. Kondo, M. Tanabashi and K. Yamawaki, Prog. Theor. Phys. 89 (1993) 1249. 

[16] S. Hands, A. Kocic, and J.B. Kogut, Ann. Phys 224 (1993) 29. 

[17] J. Schwinger, Phys. Rev. 82 (1951) 664. 

[18] I.S. Gradshtein and I.M. Ryzhik, Table of Integrals, Series and Products (Academic Press, 
Orlando, 1980). 

[19] A. Chodos, K. Everding, and D.A. Owen, Phys. Rev. D42 (1990) 2881. 

[20] P. Lee and S. Shenoy, Phys. Rev. Lett. 28 (1972) 1025. See also Z. Tesanovic, Physica C220 
(1994) 303 for an update and the list of references. 

[21] R. Haag, Nuovo Cimento 25 (1962) 287. 

[22] See the first paper in Ref.l; A.J. Niemi and G.W. Semenoff, Phys. Rev. Lett. 51 (1983) 2077; 
Phys. Rep. 135 (1986) 99; D. Boyanovsky, R. Blankenbecler and R. Yahalom, Nucl. Phys. 
B270 (1986) 483. 

[23] A. Kovner, B. Rosenstein and D. Eliezer, Nucl. Phys. B350 (1991) 235; A.Kovner and D. 
Eliezer, Int. J. Mod. Phys. A7 (1992) 2775. 

[24] N.N. Bogolyubov, Physica 26 (1960) 1; R.V. Lange, Phys. Rev. 146 (1966) 301. 

[25] D.J. Gross and A.Neveu, Phys. Rev. D10 (1974) 3225; E. Witten, Nucl. Phys. B145 (1978) 
110. 

[26] N.D. Mermin and H. Wagner, Phys. Rev. Lett. 17 (1966) 1133; S. Coleman, Commun. Math. 
Phys. 31 (1973) 259. 

[27] J. Zak, Phys. Rev. 134 (1964) A1602; J.E. Avron, I.W. Herbst and B. Simon, Ann. Phys. 114 
(1978) 431. 

[28] V.P Gusynin and V.A. Miransky, Mod. Phys. Lett. A6 (1991) 2443; Sov. Phys. JETP 74 
(1992) 216; V.A. Miransky, Int. J. Mod. Phys. A8 (1993) 135. 

[29] L. Dolan and R. Jackiw, Phys. Rev. D9 (1974) 3320. 

35 



A.M.Polyakov and P.B.Wiegmann, Phys. Lett. 131B (1983) 121. 
[31] D. Cangemi, E. D'Hoker and G.Dunne,Piys. Rev. D51 (1995) R2513. 



36 



Figure Captions 

Fig. 1. The thermodynamic potential Vp as a function of a at different temperatures at sub- 
critical g: g/g c = 0.9. All quantities are measured in fi = A/g c units: Vp — > irVp/fi 3 , eB — ► eB/fi 2 , 
a — > a J ' \x. 

Fig. 2. The thermodynamic potential V# as a function of o" at different temperatures at super- 
critical g: g/g c = 10 3 . All quantities are measured in fj, = A/g c units: Vp — > irVp/fi 3 , eB — *■ eB/fi 2 , 
a — »■ cr//i. 

Fig. 3. The critical line in the efi — T plane separating the asymmetric (A) and symmetric (S*) 
phases at subcritical g: g/g c = 0.9. All quantities are measured in /i = A/g c units: eB — > eB/fi 2 , 
T -> 7> 

Fig. 4. The critical line in the eS — T plane separating the asymmetric (A) and symmetric (S*) 
phases at supercritical g: g/g c = f.f . All quantities are measured in /i — A/g c units: e_B — > eB/fi 2 , 

t -> t> 

Fig. 5. Diagrams in next-to-leading order in 1/JV. A solid line denotes the fermion propagator 
and a dashed line denotes the propagators of er, r and n in leading order in 1/iV. 



37 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-th/9407168v2 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-th/9407168v2 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-th/9407168v2 



This figure "fig2-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-th/9407168v2 



